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Abstract: We have observed Centaurus A with the MID-infrared Interferometric in- 
strument (MIDI) at the Very Large Telescope Interferometer (VLTI) at resolutions of 
7 15 mas (at 12.5 (im) and filled gaps in the (u,v) coverage in comparison to earlier 
measurements. We are now able to describe the nuclear emission in terms of geometric 
components and derive their parameters by fitting models to the interferometric data. 
With simple geometrical models, the best fit is achieved for an elongated disk with flat in- 
tensity profile with diameter (76 ± 9)mas x (35 ± 2)mas ((1.41 ± 0.17)pc x (0.65 ± 0.03)pc) 
whose major axis is oriented at a position angle (PA) of (10.1 ± 2.2)° east of north. A 
point source contributes (47 ±11) % of the nuclear emission at 12.5 /im. There is also 
evidence that neither such a uniform nor a Gaussian disk are good fits to the data. This 
indicates that we are resolving more complicated small-scale structure in AGNs with 
MIDI, as has been seen in Seyfert galaxies previously observed with MIDI. The PA and 
inferred inclination i = 62.6^2 6 ° °f the dust emission are compared with observations of 
gas and dust at larger scales. 
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1 Introduction 

In the unified mo del of Active Galactic Nuclei (AGNs) 
( Antonucci|1993 1, the central engine (black hole + ac- 
cretion disk + broad line region clouds) is embedded 
in a dusty obscuring torus. Galaxies where we see 
the center only through this torus are called "type 2" 
whereas ones where we receive direct radiation from 
the broad line region are referred to as "type 1". For 
a long time, there was only indirect evidence for the 
existence of the enshrouding dust, until VLTI obser- 
vations with MIDI actually resolved the parsec-scale 
AGN heated dust structures in NGC 1068 ( jjaffe et al.] 
20041. Now, "tori" have been observed intcrfcromctri- 



cally in a number of Seyfert 1 and 2 galaxies ( Tristram 
et al.||2007| |Raban et al.||2009| |Tristram et al.||2009| 



Burtscher et aL][2009|) 



Radio galaxies can also be classified in type 1 (e.g. 
Broad Line Radio Galaxies, BLRGs) and type 2 (e.g. 
Narrow Line Radio Galaxies, NLRGs) sources and are 
also unified by means of orientation dependent obscu- 
ration, although a "torus" is not required in all cases 
flUrry fc Padovani|1995] ). 

Centaurus A (NGC 5128) is the nearest radio pow- 
erful AGN with a distance of only 3.80 ± 0.08 Mpc 
(see Harris et al. in this volume). At this distance, 
one parsec corresponds to 54 milli-arcseconds (mas). 



Centaurus A is therefore the only radio galaxy (clas- 
sification: NLRG) that can be studied with MIDI in 
great detail. 



With limited (u, v) coverage, Meisenheimer et al. 
( |2007[ ) found from earlier MIDI data that the mid- 
IR emission from the central parsec of Centaurus A 
is dominated (60% at at 13 /jm) by an unresolved 
point source (diameter < 10 mas). By comparison with 
multi-wavelength data, this emission was classified as 
being most likely non-thermal (synchrotron) in origin. 
The resolved emission was estimated to be > 30 mas 
(0.6 pc) and < 12 mas (0.2 pc) in size, roughly perpen- 
dicular and along the axis of the radio jet respectively 
with the major axis at about 127 ± 9 °, roughly per- 
pendicular to the radio jet. It was interpreted as a 
geometrically thin, inclined dusty disk. With only 4 
(u, v) points, the PA of the inferred disk was simply 
fixed to the direction orthogonal to the PA of the max- 
imum visibilities. 



Here we report first results from a more exten- 
sive set of (u, v) points that allows to fit model source 
brightness distributions to the nuclear mid-infrared emis- 
sion and determine the parameters of the possible emis- 
sion components more precisely. Further analyses re- 
garding the structure of the emitter will be published 
elsewhere. 
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2 Instrument, observations and 
data reduction 

Observations were perfor med in the N band (8^tm < 
A < 13/im) with MIDI (|Leinert etln~1|2003[) at the 



NGC51 28 



ESO VLTI on Cerro Paranal, Chile, using pairs of 
8.2 m Unit Telescopes (UTs). The observables with 
MIDI are the single-dish spectra and the correlated 
flux spectrum that is obtained from the interference 
pattern generated by the two beams. MIDI observa- 
tions do not provide visibility phases since the phase 
information is destroyed through atmospheric turbu- 
lence. Since MIDI can only combine the light of two 
telescopes, closure phases are also not available. 

The light was dispersed with a NaCl prism with 
spectral resolution R = A/ A A ~ 30 and a slit width 
of 200-^tm that corresponds to 0.52 arcsec on the sky. 
Since Cen A appears as a point source in the mid- 
infrared with the single-telescope observations, the ef- 
fective "field of view" of the interferometer observa- 
tions corresponds to the width of the point-spread func- 
tion at the wavelength of observation (380 mas, or 7 
pc, at 12.5 /im). 

For the observation log as well as the observational 
parameters, see Tab. [T] In total, 20 fringe track ob- 
servations were made of which one was discarded due 
to its very low signal. For another one no single-dish 
observation was taken, leading to a total of 18 visi- 
bilities. Some of them (6) are practically duplicates 
in the (u, v) plane and were taken to increase signal- 
to-noise. The resulting (u, v) coverage is displayed in 
Fig. [I] These observations provide effective spatial res- 
olution:^] at 12.5/im of 6.7 - 15.5 mas, see also Tab.[l] 
Calibrators were selected to be close in airmass. Due 
to the orientation of the VLTI baselines and the dec- 
lination of the source, there are unfortunately no long 
UT baselines in the lower left (SE) quadrant of the 
(u, v) plane. 

Data reduction was performed with the interfero- 
metric data reduction software MIDI Inter active Anal - 
ysts and Expert Work Station (MIA+EWS, |Jaffe|2004| P1 
and the data reduction procedure was as follows: First 
the associated calibrator (see Tab.[T| was reduced. For 
the extraction of the spectra, a 4-pixel wide Gaussian- 
shaped mask was used. The vertical shift of the mask 
was determined by cross-correlating the mask with the 
spectrum. Additionally, a skymask was used to de- 
termine the sky background at equal, short distances 
above and below the spectrum. 

With these masks, the science data were reduced 
with identical settings and calibrated using the cali- 
brator star. The calibrated correlated fluxes were de- 
termined directly from the raw correlated fluxes of the 
science and calibrator sources without using the single- 
dish spectra. This substantially lowers the measure- 



applying Rayleigh's criterion to an interferometer leads 
to a resolution of X/2BL. However, given sufficient S/N, 
one can already distinguish models at lower resolutions, in 
our case at ~ X/3BL. 

2 The nightly build from Feb 07 2010 

was used. It can be downloaded from 

http://www.strw.lcidcnuniv.nl/~koehler/MIA-l-EWS- 
daify-snap shot.html| 




20 -20 
U [cycles/orcsec] 

Figure 1: (u,v) coverage of the 2008 MIDI visibil- 
iti es of Cen A at 12.5 \im. Ob servations that led 
to Meisenheimer et al. (2007) (4 visibilities) are 
shown as open diamonds, the 2008 observations 
(18 visibilities) as filled diamonds. Every visibility 
appears twice in this plot. The solid lines are the 
(u, v) positions traced by the various baselines (UT 
combinations are labelled) as a result of earth's ro- 
tation. They are truncated at telescope elevations 
of 30 degrees above the horizon. 



ment error for weak sources since the uncertainty of 
the single-dish fluxes is much larger than that of the 
correlated fluxes. The intrinsic calibrator spectra were 
taken from Cohen et al.| (19991. A more detailed de- 
scription of the general MIDI data reduction procedure 
can be found in 1 



Tristram 



(20071 



All Cen A single-dish spectra were stacked and 
combined to one average spectrum (for all observa- 
tions). This is possible since only the correlated spec- 
trum depends on the projected baseline length and ori- 
entation and it is not expected that the N band source 
changes within one week. The 2005 observations that 



led to Meisenheimer et al. ( 2007 1 were not included in 



this analysis since their total flux differs considerably 
from the 2008 observations. 

Finally, the calibrated correlated spectra were di- 
vided by the average single-dish spectrum to derive vis- 
ibilities. For this work, the visibilities between 12.3/im 
and 12.7/^m were averaged and all model fits in this 
paper refer to this part of the N band spectrum. The 
full spectral information will be exploited in a subse- 
quent paper. The error of the calibrated visibility is 
composed of the statistical error of the correlated flux 
and that of the averaged photometry as well as the 
error of the calibration template (assumed to be 5 %). 
The errors of the averaged 12.5 /im visibilities are then 
on the order of 3 %, see also Tab. [I] 
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Figure 2: Visibilities in the (u, v) plane. Each sym- 
bol represents a visibility measurement, the color 
of the symbol shows the visibility amplitude. The 
dots denote radii of 20, 40 ... 140 m around (u , v) 
= 0. The data from '. 



not shown in this p 



Meisenheimer et al. ( 2007 ) are 



ot. 



3 Results and modelling 

3.1 Visibilities 

The resulting averaged visibilities at 12.5 pm (AA = 
0.2/im), measured at baselines of about 50 - 130 me- 
ters, range from about 40% - 70%. They are given in 
Tab. [I] and visualised in a (u, v) plot in Fig. [2] where 
each symbol represents a fringe observation. The color 
of the symbol denotes the visibility. Their distribu- 
tion in the (u,v) plane does not intuitively show what 
structure they probe. 

3.2 Geometric models for the sur- 
face brightness distribution 

With phase-less data on sparsely sampled (u, v) planes 
(see Fig. [I]), it is impossible to directly reconstruct 
meaningful images, i.e. images that show more than 
just the properties of the synthesized beam. We there- 
fore model the source brightness distribution (the "im- 



age") with simple geometrical components (e.g. point 
sources, rings, ellipses etc.) and constrain the param- 
eters of the model with the observed visibilities (e.g. 
Berger fc Segransan|2007 l. 



Here we present our model fit to the visibilities at 
one wavelength. The fit was performed with Lyon's 
Interferometric Tool prototype (LITpro), an interfero- 
metric model fitting tool, provided by the Jean-M arie 
Mariotti Center (JMMC, |Tallon-Bosc et al.|2008| p] It 
is based on a modified Levenberg-Marquardt algorithm 
and runs a grid-search to find the least x 2 value in the 
parameter space (in given boundaries). 



3.2.1 Two-components models 



Motivated by the finding of Meisenheimer et al. (2007) 



that the infrared nuclear emission is composed of a re- 
solved and an unresolved component, we fitted two- 
component models with a point source (unresolved) 
and an extended component (resolved). We model 
the extended source either with a circularly symmet- 
ric Gaussian intensity profile (model a), an elongated 
Gaussian intensity profile (model b) or with an elon- 
gated uniform (flat intensity profile) disk (model c). 
A circularly symmetric flat profile component, a ring- 
shaped emission component as well as limb-darkened 
profiles lead to unphysical parameters. 

All models are characterized by the flux fraction of 
the point source (/ P ). 

The size parameter has model-specific meanings: 

model a Gaussian Full- Width at Half Maximum (FWHM) 
model b minor axis FWHM 
model c minor axis diameter 

The models with an elongated extended compo- 
nent (b, c) additionally have the parameters axis ratio 
p and major axis position angle PA in degrees east of 
north. 

The center positions of the two components were 
held fixed on top of each other. Models with non- 
concentric components will be discussed in a future 
paper. 

The best fit (in terms of \ 2 ) can be achieved by 
fitting the resolved component as an elongated disk of 
(1.41 ±0.17) x (0.65 ±0.03) pc. The major axis of this 
structure is almost north-south at a PA of (10.1±2.2)°. 
In this model, the point source contributes (47 ± 11) 
% to the total flux at 12.5 /im. The parameters of 
all models are compiled in Tab. [2] Fig. [3] shows cuts 
through the \ 2 plane of the best-fitting model (c) to 
show how well constrained the parameters are: While 
the minimum for the point source flux contribution is 
well defined, the \ 2 planes show that PA and axis ra- 
tio p are not so well defined. The \ 2 plane for PA and 
p shows a number of local minima, indicating degen- 
eracies in the model fit. Fig. [4] shows a sketch of this 
model. 

3 The online version of LITpro can be found at 
|http: / /www.jmmc.fr/litpro_page.htm| 
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Table 2: Two-component model fits. For the de- 
scription of the parameters, see text. For each 
model, the % 2 value for the best fitting parame- 
ters is given together with the number of degrees 
of freedom (NDOF). 





a: Gauss 


b: El. Gauss 


c: El. disk 


fp 


0.48 ±0.11 


0.49 ±0.12 


0.47 ±0.11 


/ mas 


38.7 ±3.0 


22.9 ± 1.6 


35.2 ± 1.5 


e / pc 


0.72 ±0.06 


0.42 ± 0.03 


0.65 ±0.03 


p 




2.94 ±0.53 


2.17 ±0.17 


PA/° 




14.9 ±2.4 


10.1 ± 2.2 


NDOF 


15 


13 


13 


x 2 


135.0 


97.6 


84.1 




2 



100 




20 25 30 35 40 45 50 



e 



Figure 3: \ 2 cu ^ s 01 our model c (elongated disk). 
The color denotes the value of \ 2 (white: max- 
imum; black: minimum, see color wedge). Top: 
flux fraction of the point source versus axis ratio 
p\ Middle: major axis position angle versus axis 
ratio p\ Bottom: major axis position angle versus 
minor axis diameter O. See text for discussion. 



4 Discussion 



The parameters of ou r disk are different from the |Meisej i- 
heimer et aT] ( |2007[ ) result: They find a > 0.6pcx < 
0.2pc Gaussian disk at PA ~ 127°. The probable origin 
for this difference is their sparse (u, v) coverage. Had 
we only had two pairs of observations, e.g. sl/s2 of 
2008-04-18 a nd s9/s!0 of 2008-04-21, our results would 
be closer to Meisenheimer et al. (20071: Converting 



the visibilities of these observations to one-dimensional 
Gaussian widths, leads to full-widths at half-maxima 
of roughly 10 mas (0.2 pc) along PA ~ 30° and 20 mas 
(0.4 pc) along PA « 170°. This shows how important 
it is to study the circumnuclear dusty component of 
AGNs not only with a few baselines but with many 
different telescope combinations. 

In fact, the visibilities look rather chaotic with val- 
ues as high as 70% separated by only about two tele- 
scope diameters from values as low as 55%. Possi- 
ble causes for these variations, other than the effects 
from the science source itself, could be variations of the 
atmospheric transfer function between calibrator and 
science target. Typical variations of the transfer func- 
tion were « 3% in between science and calibrator ob- 
servation as determined by linear interpolation of the 
transfer function between each two subsequent calibra- 
tor observations. Another possibility for artificial flux 
variations are differences in the beam profile between 
observations (e.g. due to vignetting). This effect could 
be important since we did not divide each correlated 
flux spectrum by its corresponding single-dish spec- 
trum individually but averaged all single-dish spectra 
in order to reduce noise (see Section|2|. We tested this 
by looking for differences in the calibrator single-dish 
beam profiles and found none of significance. Further- 
more, it is possible that spectra observed on different 
baselines have systematic offsets. We tested this by 
observing the same (u,v) point twice on two different 
baselines. The scatter between the baselines was not 
larger than the scatter within the baseline. 

These possible causes for artificial visibility varia- 
tions cannot explain the differences observed. There- 
fore we believe these offsets are real. 
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VLTI resolution 
(~ 7 mas) 




taken into account that could artificially create elon- 
gated disks. This can possibly be constrained with 
future observations. 

The highest effective spatial resolutions of the new 
12.5 fim data presented here are still about 10 x larger 
than the sublimation radius of dust in Cen A ( jMeisen-] 
heimer et al. 12007] ) and the point source of our model 
could therefore still be contaminated by dust emission 
from the innermost part of the dusty disk. A detailed 
spectral model, taking into account also the visibili- 
ties at other wavelengths, will be able to decompose 
the two components more precisely. 



Figure 4: Two-component model for the nuclear 
mid-infrared emission of Centaurus A consisting of 
a uniform disk component and an unresolved point 
source that is smaller than the VLTI resolution 
(dark grey circle) and has been identified with the 
synchrotron core (black point) by Meisenheimer 
et al.| ( |2007"| . The resolved structure is ~ 1.4x 0.65 
pc at 12.5 /jm and the best fit PA ~ 10.1°; the 
radio radio jet axis is PA ~ 51° (e.g. Tingay et al. 



2001) 



In other sources such jumps in the visibility are 
only seen at very low visibility levels (< 20 %) where 
the source is almost c ompletely resolved (cf. Fig. 4 of 
Tristram et al. ( 2007 1). But in Cen A the point source 
contributes almost 50% to the total flux according to 
our model fits. Therefore, we are effectively observing 
the extended component at visibilities < 20 %. This 
means we are probing the extended component at rela- 
tively high spatial resolutions which are most sensitive 
to small-scale structures and not well described by a 
Gaussian or uniform disk. 

This is not only seen in our Cen A data, where 
our best fit model results in a \ 2 — 84.1 (with 13 de- 
grees of freedom), but also in the two other AGNs that 
have been studied extensively with MIDI: The Circi- 
nus galaxy and NGC 1068. In the Circinus galaxy, a fit 
with two Gaussian compone nts leads to a \ 2 ~ 1 6600 
with 451 degrees of freedom (Tristram et al. 1120071); in 



NGC 1068 especially the data at longer baselines are 
also not well described by a mo del of two Gaussian 
components (Raban et al. 20091. These small scale 



structures could, for example, arise from dumpiness 
in the dusty disk as has been demonstrated by |Tris-| 
|tram et al.| ( |2007| ). 

It is therefore not completely unexpected that the 
geometric fit to the Cen A data doesn't fully describe 
the visibility data. 

A comparison of our models a, b, c shows that (i) 
an elongated disk considerably improves the fit and 
(ii) a Gaussian profile leads to a similar disk area (~ 
67mas x 23mas, corresponding to ~ 1.2 pc x 0.4 pc) as 
a flat profile (~ 76mas x 35mas, corresponding to ~ 1.4 
x 0.6 pc) at very similar position angles of (14.9±2.4)° 
and (10.1 ± 2.2)°, respectively. However, when inter- 
preting the reliability of the PA of the disk in such fits, 
the asymmetric (u, v) coverage (Figs, [I] |2| must be 



4.1 Comparison with other observa- 
tions 



Neumayer et al. (20071 used adaptive optics assisted 



integral field spectroscopy in the near infrared to derive 
a model for the mass distribution of the central 3 arcsec 
(60 pc) from the observation of molecular hydrogen 
emission. At their resolution limit of 120 mas (ca. 2 
pc) in the K band, they find a PA of (148.5 ± 1)° and 
an inclination of (37.5 ± 2)° for the velocity field in 
their tilted ring model. 



Espada et al. ( 2009 \ observed the central 1 arcmin 
(1 kpc) of Cen A in the 12 CO(2 — 1) emission line with 
the Submillimeter Array with a resolution of 100 pc x 
40 pc. Using a warped thin disk model, they find that 
this molecular gas emission in the inner 20 arcsec (400 
pc) is elongated at a PA of ~ 155° and they interpret 
the observed axis-ratio of the presumably circular disk 
as an inclination of i — 70°. In VLBI observations, 



Tingay et al. ( 2001 1 found a position angle of the jet of 
51 and confirmed the previously found inclination of 
50° - 80° ( |Tingay et al.|1998|). In combined radio and 



X-ray observations, Hardcastlc ct al. (2003) find that 



the radio and X-ray jet trace each other quite closely 
in terms of position angle and that only small incli- 
nations (angles to the line of sight) of « 20° are con- 
sistent with both X-rays and radio observations from 
combined constraints on apparent motions and sided- 
ness of the jet (components). The PA of the X-ray jet 
is ~ 55° (Kraft et al.|2000b. 



Assuming that the dusty component is well rep- 
resented by an inclined disk, our axis ratio of p — 
2.17 ± 0.17 can be interpreted as an inclination of 
cos -1 1/p = 62.6I2 6 ° an d w ould then lie roughly in 
bet ween the measuremen ts of Neumayer et al. ( 2007 1 
and |Espada et al.| ( [2009 ). 

The position angle of ~ 10° (= 190°) from our fit 
is rather unexpected when compared with these ob- 
servations on larger scales, however. From the naive 
unification picture argument, one might expect the gas 
and dust to be elongated along ~ 140°, at a right an- 
gle to the radio and X-ray jet. This differs by ~ 50° 
from our model fit. Before explaining this with physi- 
cal effects (warped accretion disk?), however, we need 
to better understand the chaotic visibility pattern in 
our observations to better constrain the position angle 
of the resolved component. 
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5 Summary 

We presented new interferometric data from MIDI/VLTI 
observations that allow us to constrain the geometry 
of the mid-infrared emission at 12.5 u,m from the nu- 
cleus of Centaurus A. We find that, in the context of 
simple geometrical models, this emission is best re- 
produced by a two-component model of an unresolved 
point-source (< 7 mas) and an elongated (Gaussian or 
fiat) disk with a size of about 1.3 x 0.5 pc (average 
of models b and c). Assuming that the dusty compo- 
nent is well represented by an inclined disk, the disk 
has a PA of about 12.5° and is inclined by about 66° 
(averages of models b and c) 

The new data have also shown that the extended 
component is almost completely resolved and that we 
are probing it at very low visibilities / high spatial 
resolutions - and that we probably resolve small scale 
structure that is not described by our geometrical model. 
The fit parameters for the disk should therefore be in- 
terpreted only as tentative evidence for such a struc- 
ture. 
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